Interference effect on Raman spectrum of graphene on Si02/Si 
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Abstract 

The intensity ratio between two major Raman bands in graphene is one of the most important 
information for physics of graphene and has been beUeved to represent various intrinsic 
properties of graphene without critical assessment of extrinsic effects. We report a micro 
Raman spectroscopy study on the Raman intensity ratio of the 2D band to the G Raman band 
of graphene varying the thickness of dielectric layers (SiOa) underneath it. The ratio is shown 
to change by almost 370% when the thickness is varied by 60%. The large variation in the 
ratio is well explained by theoretical calculations considering multiple Raman scattering 
events at the interfaces. Our analysis shows that the interference effect is critical in extracting 
the intrinsic 2D to G intensity ratio and therefore must be taken into account in extracting 
various physical properties of graphene from Raman measurements. 
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I. INTRODUCTION 

Graphene, a two-dimensional hexagonal crystal of carbon atoms, has attracted immense 
interests from researchers in various disciplines because of its novel electronic properties, 
such as a high carrier mobility^"^ and anomalous quantum Hall effect."*'^ These intriguing 
properties are caused by the linear energy dispersion versus momentum around the Dirac 
points/"^ After the first successful isolation of graphene,^ its unique physical properties have 
been studied by using various experimental tools: two or four terminal transport 
measurement,"^'^ Raman spectroscopy,^ spin transport,^ infrared spectroscopy,^ angle resolved 
photoemission spectroscopy,^ and scanning tunneling microscopy/ 

Among these experimental probes, Raman spectroscopy is one of the most successful tools 
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in investigating the electronic and structural properties of graphene. ' " A typical Raman 
spectrum of graphene consists of two major features (Fig. 1): the G (cog ~ 1586 cm"') and the 
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2D ((02D ~ 2686 cm" ) bands. ' ' The G band originates from the Stokes Raman scattering 
with one phonon {E2g) emission. It is known that as the doping concentration is increased, 
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its frequency blueshifts and its width decreases. " The 2D band (or sometimes called G') is 
due to the Stokes-Stokes double resonant Raman scattering with two-phonon (Ai') 
emissions,^''^ and its shape is very sensitive to the number of graphene layers. ^''^"''* Although 
the absolute intensity of either of these peaks in a Raman spectra may depend on various 
external factors such as equipment alignment that may vary in each measurement, the 
intensity ratio of the 2D band to the G band (I2d/Ig) is often thought to be immune to such 
external factors and represent the intrinsic properties of a given type of graphene. It has been 
used to determine basic structural and electronic properties of graphene such as the number of 
layers,^''^"'"* doping concentration,'*'"^' and optical anisotropy.^^ 
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In this Letter, we find that through micro-Raman spectroscopy, ' the observed ratio 

Iil/^g varies by 370 % when the thickness of the SiOi layer on silicon substrates is varied by 

2 



60 %. Hence, the thickness of the SiOi layer, which should not affect the intrinsic properties 
of the graphene sample on it, should be taken into account in interpreting the observed hr/Io- 
By considering the interference of the excitation laser as well as the Raman signal due to 
multiple reflections at the interfaces, one can calculate the enhancement factors for the 2D 
and the G bands. We show that these factors are different in general, and therefore, one needs 
to factor out these enhancement factors properly when important intrinsic properties such as 
Iid/Ig are deduced from the experimental data. 

II. EXPERIMENTAL 

Mechanically exfoliated graphene samples^ were placed on top of SiOa/Si substrates with 
various SiOi layer thicknesses (~240 to ~380 nm). Initially, Si (p-type) substrates covered 
with ~300- or ~388-nm SiOa layer were prepared by wet thermal oxidation. Then, the 
thickness of the SiOa layer was reduced by wet etching in a buffered NH4F-HF (BHF) 
solution for various etching times. The surface roughness of unetched Si02 is about 0.16 nm 
as measured by atomic force microcopy. Etching increases the roughness, but all the etched 
substrates have more or less similar roughness (0.54—0.62 nm). 

Since the thickness and refractive index of the Si02 layer are crucial factors in the analysis 
of the enhancement factors for the Raman intensity, we used high-precision spectroscopic 
ellipsometry (SE) to determine their precise values. The SE measurement^"*'^^ was performed 
under high purity N2 atmosphere in the wavelength range from 190 to 1100 nm. To reduce 
experimental errors we measured the SE spectra at multiple incidence angles of 60, 70, and 
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80°, and then exfracted a single set of parameters to fit all 3 spectra. The error bars in the 
determination of the thickness and the refractive index are ±0.3-2.0 nm and 2.0x10"^, 
respectively. 
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Single layer graphene samples are roughly identified with an optical microscope, and then 
confirmed by micro-Raman spectroscopy measurements. Because of the interference effect 
between the graphene and SiOa layers, the color and the contrast of graphene are influenced 
by the wavelength of the illumination and the thickness of the SiOi layer. " In general, it is 
known that a ~300-nm thick SiOi layer is optimal for white light illumination. Since 
graphene samples on SiOa layers thicker than 340 nm are not visible under white light 
illumination, a red dichroic filter was used for these samples (Fig. 1 (a) and (b)). 

For the micro Raman measurements, the 514.5-nm (2.41 eV) line of an Ar ion laser was 
used as the excitation source, and the laser power was ~1 mW. The laser beam was focused 
onto the graphene sample by a 40x microscope objective lens (0.6 N.A.), and the scattered 
light was collected and coUimated by the same objective. The scattered signal was dispersed 
with a Jobin-Yvon Triax 550 spectrometer (1200 grooves/mm) and detected with a liquid- 
nitrogen-cooled CCD detector. The spatial resolution was less than 1 |xm, and the spectral 
resolution was about 1 cm"\ The single layer graphene was exactly identified by the unique 
shape of the Raman ID band as shown in Fig. 1 



III. RESULTS AND DISCUSSION 

We found that the Raman intensities of the G and 2D bands are indeed strong functions of 
the thickness of SiOa layer [Fig. 2(a)]. As the thickness is increased from 240 nm, the 
observed intensities increase first and show the highest values at the thickness of ~280 nm for 
the G band and ~290 nm for the 2D bands, respectively. Since the two maxima occurs at 
different wavelengths, the resulting observed ho/Io ratio varies greatly; the maximum of the 
intensity ratio is around 9.3 and 6 times higher than the minimum (Fig. 2(b)). Some scatter in 
the experimental data will be discussed later. 



4 



To explain the observed peculiar variation of the intensity, we use the multi-reflection 
model (MRM) of the Raman scattered light (Fig. 3 (a) and (b)). In this model, the absorption 
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and scattering processes are treated separately. We note that a similar method was applied 
to explain the variation of the G band intensity when the number of graphene layers 
increases. We also note that it was applied schematically to Raman intensity variation as a 
function of the thickness of the dielectric layer without considering the difference in the 
wavelengths of the laser and Raman scattered light. In general, the wavelengths of the 

3 1 

Stokes Raman scattered light and the laser are taken to be the same in similar calculations. 

However, in graphene, the actual differences between the wavelengths of the laser, the 
Raman G band, and the 2D band are quite large. When the 514.5-nm (2.41 eV) line of an Ar 
ion laser is used, the Raman G and 2D bands of single layer graphene are located at ~1586 
cm"^ and ~2686 cm"\ respectively. In terms of wavelengths, these correspond to 560.2 nm 
(2.21 eV) and 597.0 nm (2.08 eV), respectively. Hence, these differences and concomitant 
differences in the index of refraction will result in different interference patterns for the 
Raman bands. 

In the MRM, as shown in Fig. 3(a), the laser beam is absorbed by the 7r-electrons of 
graphene while passing through the graphene layer. However, the laser beam goes through 
multiple reflections inside the graphene layer as well as in the Si02 layer. Due to these 
multiple reflections, there are multiple chances for the beam to be absorbed by the k- 
electrons. The net absorption term (Fab) could be represented by the sum of the dots in Fig. 
3(a) and can be expressed as 
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where = In^ / (n^+n^) , r^= («(, - «,) / («q + «,) , r^={n^ — n^) I (n^ + "2) 

=(^2 -«3)/(«2 + "3) the Fresnel transmittance and reflection coefficients for the 
interfaces involving air (0), graphene (1), Si02(2), and Si (3). is the refractive index 

of air, and n^ , and are the refractive indices for graphene, SiOi, and Si, 
respectively. We also use abbreviations fi^ = 2n:xn^ / A, Ind^h^ IX, P2- '^^d^n^ I X , 
where x is the depth of the point where the interaction occurs, and and d.^ are the 
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thickness of the single layer graphene and the SiOa layer, respectively. 

Similarly, the net scattering term (F^c) could be represented by the sum of the arrow lines 
in Fig. 3(b) and expressed as 



where t[ — 2n^ /(n^+n^). A is the wavelength of the excitation source in the net absorption 
term and is the wavelength of the G or 2D bands in the net scattering term. Then, the total 
enhancement factor (F) is given by 



where TV is a normahzation factor, which is a reciprocal number of the total enhancement 
factor for a free-standing graphene, obtained by replacing the SiOa and Si layers with air. The 
measured Raman intensity (I) is I=Ii-F, where /, is the intrinsic Raman intensity playing the 
role of a single fitting parameter in the subsequent calculations. 

Figure 3(c) is the calculated enhancement factor (F), relative to the case of a free-standing 
graphene film. The thickness of graphene ( Jj ) is taken to be 0.34 nm, which corresponds to 

the interlayer distance in graphite crystals. The interference effects on the Raman G and 2D 
bands are clearly seen; as a function of the thickness of the SiOa layer, the enhancement 





(3) 
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factor for the Raman G and 2D bands vary by a factor of up to 48. In Fig. 2(a) and (b), the 
results from our model calculation are compared with the experimental data for the Ig, ho, 
and hr/Io- Overall, they show good agreement with each other. We can fit the data in Fig. 
2(a) by setting the intrinsic Raman intensity ratio //,2Z)/ h.G =3.4. Although there is some 
scatter, the measured data are well represented by the curves fitted from the model 
calculation. It should be noted that the calculated enhancement without considering the 
difference in the wavelengths for each Raman scattering event (solid curve in Fig. 3(c)) 
deviates significantly from the correct ones when the dielectric layer becomes thicker than 
200 nm. 

We found that some of the scatter in the data for the intensity ratio originate from other 
extrinsic factors such as defects and doping. '^'^^ Some of the samples showed the defect- 
induced Raman D band (-1350 cm"^) signals. These samples tend to have higher G band 
intensities as compared with those that do not show the D band. In Fig. 2(b), the data from 
these samples (star symbols) are off the main tendency. Other samples showed higher levels 
of doping, as indicated by blueshifted G band peak positions and decreased widths of the G 
band."^"^^ Highly doped samples are known to exhibit lower Ixl/Ig ratios. '^'^^ The data from 
the samples with estimated doping densities'^"^^'^^ in excess of 5x10'^ cm"^ are identified 
by open circles in Fig. 2(b). By excluding data from such samples, we fit the data from 
remaining 'intrinsic' samples to obtain /,-2z)/ Ag=3.2. 

In the above analysis, we used a value of 0.34 nm as the thickness of single layer 
graphene. However, the thickness of graphene, which is just one atomic-layer thick, is not a 
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well-defined quantity. AFM measurements ' do not give an definitive answer as to the exact 
value of the 'thickness' of graphene. Also it is known that there exist ripples on the order of 
about 0.5 nm, which also affect the 'optical thickness' of the graphene layer. In order to test 
the sensitivity of our model to the choice of the thickness value, we repeated the calculation 
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while varying the thickness of single layer graphene from 0.1 to 1 nm. However, this changed 
the interference pattern only slightly: the maxima of the enhancement factors (Fig. 3(c)) 
shifted by less than 2 nm. There are also speculations that there exists a thin layer of air or 
water between graphene and Si02, which might also affect the interference. Again, we 
calculated the interference pattern as a function of the thickness of interlayer (air or water) in 
the range of 0-1 nm, which shifted the maxima of the enhancement factors by less than 2 nm. 
The index of refraction of graphene is another ambiguous quantity. We used the index of 
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refraction of graphite (n^ ) as a first approximation as other authors did. " This somewhat 
arbitrary choice does not affect our result, though, because varying the index of refraction 
from 0.5«j to 2.0«j changed only the absolute amplitude of the interference pattern but did 

not shift the interference pattern. As a matter of fact, the interference pattern changed 
appreciably only when the thickness or the refractive index of Si02 was varied. 

In general, the interference effect is a fiinction of the incident angle of the light. In our 
analysis, we assumed normal incidence because most micro Raman scattering measurements 
are performed in backscattering geometry. Even when the numerical aperture (N.A.) of the 
objective lens is large, the laser beam is almost normally incident on the sample provided that 
the beam is Gaussian and the focused laser beam hits the sample surface at the beam waist. 
However, in practice, since the focal depth is about 1 |j,m for N.A.=0.6 and A = 500 nm, it is 
conceivable that a significant portion of the beam enter the sample at an oblique angle if the 
focus is only slightly off. We considered the upper bound of the effect of the large N.A. on 
the interference pattern, regarding a Gaussian disfribution of the incident light intensity but 
treating the beam path with classical ray optics. Contributions from each portion of the beam 
with an incident angle (O<O<0j^ =arcsinN.A.) were calculated separately and then 
integrated over 0. As shown in Fig 4(a) and (b), the patterns of the enhancement factors of 
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the G and 2D bands shift slightly for larger N.A. values. The peak positions for 0.9 N.A. is 
shifted by about 10 nm with respect to those for normal incidence. The ratio of the 
enhancement factors {F2d^Fg) is also slightly shifted in Fig 4(c). New fitting curves including 
the effect of a large N.A. value of 0.6 are shown in Fig. 2 as dashed curves. These new fitting 
curves give Ii2d/ kcn'^.S. 

We also calculated the enhancement factors of the G (Fg) and 2D (F2d) bands as functions 
of the thickness of SiOi layer and the wavelength of the excitation source for normal 
incidence. The results are shown in Fig. 5(a) and (b). In this calculation, co^. is fixed since it 
does not vary with the laser wavelength. However, the frequency of the Raman 2D band ( 
(O2D ) depends on the laser wavelength, and the dispersion of the Raman 2D band can be 
given as a linear function, co^d - 2444.24 + 99.06£'^^^^cm"' , where E laser is the laser energy 

38 

in eV. Figures 5(a) and (b) clearly show that the Raman signal is significantly enhanced or 
suppressed, depending on the laser wavelength and the Si02 layer thickness. Figure 5(c) is a 
contour plot of the ratio of the enhancement factors for the Raman 2D band to the G band. It 
is clear that these factors play a major role in determining the intensity ratio. Our calculation 
indicate that it is important to factor out the interference effect first, when comparing ho/^G 
data fi-om samples with different SiOa layer thicknesses, or obtained with different lasers. It 
should also be noted that this kind of interference effect is not unique to Raman 
measurements but applies to any spectroscopic measurements on thin samples on dielectric 
layers. Appropriate choice of the dielectric layer thickness therefore can significantly enhance 
the measured signal in such cases. 



IV. SUMMARY 

In summary, a strong dependence of the Raman spectrum of single layer graphene on the 
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thickness of the SiOi layer on the substrate is observed and anal5^ed in terms of multiple 
reflection interference. It is found that the Raman spectrum depends not only on the Si02 
layer thickness but also on the wavelength of the excitation laser. This effect significantly 
influences the observed intensity ratio of the Raman 2D band to the G band. 
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CAPTIONS 



FIG. 1. (Color online) Raman spectrum of a single layer graphene sample taken with a 514.5- 
nm laser as the excitation source. Insets (a) and (b) are the optical microscope images of a 
single layer graphene sample on a 350-nm SiOa layer with white and red (filtered, 615-730 
nm) light illumination, respectively. The graphene layer can be easily seen in (b), but not in 
(a). 

FIG. 2. (Color online) (a) G (circle dots) and 2D (square dots) band Raman intensities as 
functions of the thickness of the SiOa layer, (b) Raman intensity ratio hr/Ie as a function of 
the thickness of the SiOi layer. The stars represent data taken fi-om Raman spectra which 
showed the D band (~1350 cm"'), and the open circles represent data taken fi-om the samples 
with high doping. The curves in (a) and (b) are the calculation results based on the MRM 
model. The inset is the calculated result for to 500 nm. The dashed curves are the results 
when the effect of the large N.A. is included for N.A.=0.6. 

FIG. 3. (Color online) Schematic diagrams of multiple reflection interference in the (a) 
absorption and (b) scattering processes. , n^, and «3 are the refiactive indices of air, 
graphene, SiOa, and Si, respectively, and d^^Q^ are the thickness of graphene and Si02 

layer, respectively, and x is the depth in the graphene layer. The dots are the points of 
interaction between the laser beam and the 7r-electrons of graphene. (c) Calculated Raman 
intensities as a fimction of the thickness of the Si02 layer. The solid curve was obtained in a 
simple approximation where the Raman bands and the laser are taken to have the same 
wavelength. 
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FIG. 4. (Color online) Calculated Raman enhancement factors of (a) G band and (b) 2D band 
as functions of the thickness of the SiOa layer for various values of the numerical aperture of 
the objective lens. The excitation wavelength is taken to be 514.5 nm. (c) Ratio of the 
enhancement factor for the 2D band to that of the G band, FjjJFq. 

FIG. 5. (Color online) Plots of calculated Raman enhancement factors of (a) G band and (b) 

2D band as functions of the thickness of SiOa layer and the wavelength of the laser, (c) Ratio 
of the enhancement factor for the 2D band to that of the G band, F^lJFg. 
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